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ABSTRACT 

The  radiation  impodanoe  of  a  finite  array  of  finite  cylinders  is 
determined  by  a  boundary  value  problem  technique.  The  final  equations 
for  the  radiation  impodanoe  are  stated  without  rostriotion  on  oylinder 
length  or  radius,  cylinder  spacing,  or  on  the  amplitude  and  phase  of 
the  radial  velocity  of  the  cylinders.  The  derivation  is  made  for  the 
case  of  each  oylinder  in  the  array  having  a  uniform  radial  velocity  over 
a  sector  of  the  oylinder,  but  the  method  of  solution  is  not  limited  to 
suoh  a  velooity  distribution.  The  results  of  a  numorical  oaloulation  of 
the  effeot  of  electrical  beam  steering  on  the  radiation  impedance  of  a 
speoific  array  are  shown  to  illustrate  such  effeot  and  the  use  of  the 
derived  equation  for  radiation  impodanoe. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance  and 
enoouragoment  given  him  by  LIr.  John  Hickman  of  the  U.  S.  llavy  Electronics 
Laboratory  and  Professor  II.  Lledwin  of  the  U.  S,  Naval  Postgraduate  School 
in  this  investigation. 
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The  notation  associated  with  the  Bessel  Functions  used  in  this  paper 
is  taken  from  N.  W.  McLaohlan,  Bessel  Functions  for  Engineers,  Oxford 
at  the  Clarendon  Press,  1955. 

Jn(z)       Bessel  function  of  the  first  kind 

Yn(z)       Bessel  funotion  of  the  second  kind,  as  defined  by  Weber 

K^z)      Modified  Bessel  funotion  of  the  second  kind 

(l),(2) 

IL^         Bessel  function  of  the  third  kind* 

H^z)  =  J  (z)  /  jY  (z) 
n       an 

(21  , 
Hn  (z)  =  Jn(z)  -  jYn(z) 

N  The  positive   square  root   of  a  real,   positive   quantity. 


1*   Introduction 

The  radiation  impedance  of  a  finite  array  of  finite  cylinders  in  an 
infinitely  long,  rigid,  cylindrical  baffle  has  been  investigated  by  D.  H. 
Robey.   Robey  considered  the  problem  of  uniformly  pulsating  cylinders 
and  his  final  conclusions  are  based  on  a  cylindrical  array  of  radius 
negligible  in  comparison  to  the  wavelength  of  the  radiated  sound.   It  is 
the  intent  of  this  paper  to  investigate  the  radiation  impedanoe  of  such  a 
finite  array  of  finite  cylinders  with  fewer  restrictions  imposed  on  the 
radial  velocity  distribution  and  size  of  the  array.   It  is  also  the  intent 
of  this  paper  to  apply  the  results  of  such  an  investigation  to  the  study 
of  the  effects  of  electrical  beam  steering  on  the  radiation  impedance  of 
a  finite  array  of  finite  cylinders. 

The  radial  velocity  distribution  on  each  pulsating  cylinder  will  be 
specified  as  a  boundary  oondition  to  the  wave  equation  which  is  to  be 
solved.  However  this  is  not  enough  to  completely  specify  a  unique  solu- 
tion to  the  wave  equation.  Therefore  it  is  proposed  that  as  the  radial 
distance  from  the  cylindrical  radiator  inoreases,  vanishing  boundary 
values  for  pressure  and  velocity  are  required.  This  condition  again 
proves  too  weak  to  determine  a  unique  solution.  Therefore  a  more  stringent 
boundary  condition  in  addition  to  the  condition  of  vanishing  boundary 
values  is  required,  ie.  "The  Sommerfeld  Condition  of  Radiation."  That  is, 
"The  energy  which  is  radiated  from  the  sources  must  scatter  to  infinity; 


D#  H.  Robey,  On  the  Radiation  Impedanoe  of  an  Array  of  Finite  Cylinders, 
J.  Acoust.  Soo.  Am. ,27,  pp706-711,  1955 

^A.  Sommerfeld,  Partial  Differential  Equations  in  Physics,  Academic 
Press  Inc.,  1949,  p  189 


no  energy  may  be  radiated  from  infinity  into  tho  prescribed  singularities 
of  the  field  (plane  waves  are  excluded  since  for  them  even  the  condition 
u  -  0  fails  to  hold  at  infinity)." 

The  wave  equation  which  is  to  be  satisfied  is  stated  in  cylindrical 
coordinates : 

6r     r   9r    r   "d  e     0  z  c   d  t 

The  scalar  velocity  potential  function  is  understood  to  be  ^(r,e,z,t). 
The  z  axis  of  the  cylindrical  system  corresponds  to  the  longitudinal  axis 
of  the  cylindrical  array  (see  figure  l). 

The  driving  forces  which  are  to  be  imposed  on  the  elements  of  the 
cylindrical  array  will  vary  harmonically  with  time.  Therefore  for  the 
results  of  this  paper  to  conform  to  standard  electrical  engineering 
practice  of  assigning  positive  reactance  to  inductance,  a  time  dependence 
of  the  form  eJ^  will  be  assumed.  That  is  jZf(r,  e,z,t)  =  0(r,e,z)  eJ**. 
After  substitution  for  0(r,ez,t),  the  wave  equation  is  reduced  to  the 
following: 

Sr2    r     dr    r2   9e2   9z2     c2 
A  formal  solution  to  the  above  equation  is  discussed  in  Appendix  1. 

The  radiation  impedance  of  an  array  of  finite  cylinders  will  now  be 
determined. 


2o   Derivation 

The  first  array  to  be  investigated  is  the  least  restricted  in 
boundary  conditions  of  the  arrays  to  be  considered,.  The  array  is  to  be 
constructed  as  follows  (see  figure  l): 

1,  There  are  m  active  cylinders  of  length  =  2^*,  radius  =  a, 

and  center  to  center  spacing  =  2b, 

2.  The  cylinders  are  mounted  coaxially  in  an  infinitely  long, 
rigid,  cylindrical  baffle  of  radius  =  a. 

The  pulsation  of  any  active  cylinder  is  to  be  described  as  follows: 
1«   On  the  given  cylinder,  the  radial  velocity  is  constant 
over  the  length  of  the  cylinder  for  a  given  aximuth  angle  a 
2.   On  all  cylinders,  the  radial  velocity  is  uniform  between  azimuth 
coordinates  -  <X   to  -f-  <K    and  zero  elsewhere. 

Therefore  the  radial  velocity  of  the  i"^^1  cylinder  may  be  described  by  the 

following  Fourier  Series: 

oo 


v  o  V,  eJ^1   /  A  cos  ne 
r  n  (EoOl) 


n  =  0 


Where     An  -  sn  sin  noc  ,     cn     a  1,        n  «  0 

ncc 

€^     ~   c i        n  ~   -L,£,o,«»»« 

and  Vi   -  the  phase  of  the  pulsations* 

The  scalar  velocity  potential  0(r,e,z,t)  must  satisfy  the  following 
wave  equation: 


2 
V2^=^    °     t  (2.02) 

at  * 


The'scalar  velocity  potential  is  assumed  to  vanish  as  r-*-oo  and  the 
Sommerfeld  Radiation  Condition  is  assumed  to  be  valid  with  respeot  to  the 
radial  coordinate. 

The  time   dependence   of  the    scalar  velocity  potential   is  taken  to  be 
j?f(r,e,z,t)  -  j^(r,e,  z)&Jwr^   as   discussed   in  Section   I.     The   scalar  velocity 
potential  ^(r,e,z)  must  then  satisfy  the  Helmholtz  V/avo  Equation 
(V  -/lc    )  0  "   0»      (2.02)     An  eigenfunotion  which   is   a  solution  to  the 
Helmholtz  Wave  Equation,    satisfies  the  vanishing  boundary  conditions   and 
the   Sommerfeld  Radiation  Condition,    is  developed  in  Appendix  1, 

This  eigenfunotion  is    stated: 

d  (  \-n  sin(ne)        -j2ir  f  z        (2).         , 

Ky    (r.e.O  -  ^   cos(ne)  lln        ^r>  (2.03) 

Y/here  /c 2  *  k2  -    (2D-  /  )2  (2#04) 


y±   =      <k2  -    (2^J  )2  (2^r)2<k2 

x*  =  -j  ^2^  f  )2  -  k2"  (2  /r'  J^  )2  >  k2 

It   is  now  proposed  to  determine  a  $  (r,e,z),   constructed  from  a  sum 
of  the  eigenfunctions  fin  v>    (r,e,z),  which  will   satisfy  the  boundary  condi- 
tions  at  the   surface   of  the  cylindrical   array.     The  technique  to  be 
employed  in  determining  the  complete   scalar  velocity  function  is  a 

■z 

logical  development  from  tho  technique  of  L.  V,  King  used  in  the  solution 

of  the  radiation  impedance  of  a  piston  in  an  infinite  rigid  baffle.  The 

constant  An  w>  transforms  into  fn(  J  )  as  the  eigenvalue  J    becomes  a 

continuous  variable. 

The  sum  of  eigenfunctions  is  taken  to  be: 
0°° 

00  \  .     . 

^(r,e,z)   «y        \  f^Cf)  cos(ne)  Hn(2j(/r)   e^T  «  df  {^Q5) 


n«0  I 

w 

-00 

3 
L.   V.  Kihg,    On  the  Acoustic  Radiation  Field  of  the     Piejfd-Electric 

Oscillator   and  the  Effect   of  Viscosity  on  Transmission,   Can.Journ,   Res.,    11, 

pp   135-155,    July-Dec. 1934 
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The  velocity  in  the  radial  direction  is: 
(r,ez)  =  -  -It-  $   (r,«,z) 


v 


oo 


vr(r,e,z)  a 


oo 


n=0 


(2) 


(2) 


fn(f)   Wl^r)  -Hn-l(/r)  cos(ne)  ej^; 


d]> 


-oo 


(2.06) 
The  radial  velocity  distribution  of  the  finite  array  will  now  be 
investigated.  From  the  boundary  conditions  as  defined,  it  is  possible  to 
write  an  expression  for  the  radial  velocity  amplitude  at  the  surface  of  the 
array: 


v  (a,e,z) 


oo 


n=0 


g(z)  An  cos(ne) 


(2.07) 

The  function  g(z)  is  equal  to  zero  on  the  rigid  baffle  and  on  the  i^*1 

cylinder  is  equal  to  V±   e^'i.   If  the  two  expressions  for  the  radial 

velocity,  (2.06)  and  (2.07),  are  compared;  the  equality  of  the  two 

quantities  is  possible  if  they  are  equal  term  by  term,  that  is: 
oo         /g)         (2) 

sU)An  .  Wj)  Vi(Ar)-^iUr)^  •»»/  -  dJ. 

-°°  (2.08) 


Define  a  new  quantity,  (3)  (2) 

'Wl^a)  -  Hn-1    (Aa) 
2 


»<r>.-r-'»(Jr> 


An 


A 


(2.09) 


Then 


•»oo 


S(z) 


G(f  )  e   4fr*V  z     dj 


(2.10) 


-CD 


It  can  now  be   seen  that   G(f  )   and  g(z)   are   a  Fourier  Transform  Pair. 
Therefore : 


..oo 


G(f)»      \  g(z)  e-J2?r^  z  dz 


(2.11) 


-oo 


Vj  oJ^l  e^2^  z  dz     -f 


-oo 


^b-£  rb//3 

(0)  e-J'2^  z  dz     y^ 
°  b-> 

(21-Db^ 

^ ^  ^  Vi  e^i  e^'2*!?  z     dz     / 

(2i-l)b~yeT 

r   (2m-l)b-//3  roe 

7^ 7^  Uej^e  -^i     d^  \        (0)e"J2^zd; 

(2m-l)b-/s  (2m-l)b^ 

sin  27r^ 


c(y)  =  v1  eJ^VJzn-y  b 


iry 


■/ 


^  ....-/  Vm  eJ^  e"J2T7t  (2n-Db         sin  2^ 


m 


a(J)  = 


i=l 


V.   ej^i  e-J27r!f  (21-l)b 


sxn 


zjrtf> 


fr 


(2.12) 


It   is  now  established,    after   substitution: 

m 

2^  »in(2rtf)    ^>   Vi  e^i  e'^^f  (2i-l)b 


In  CD  = 


i=l 


TT2l  (2l         7 


(2.13) 


and 


oo 


oo 


(2) 


/(r,e,z)= 


n=0 

L-oo 


2An  cos(ne)   sin(2^f/3  )Hn        Ur) 
TT2)  (7)         IT 


X 


m 


"\    Viejn  e  -J27r^    (2i-l)b     9j2/T^  z  d^ 


l«a- 


(2.14) 


To  obtain  the   average   pressure,    one  uses  the   foroe  equation 

vp/e^  :  ° 

Sinoe   V  ■  -  V0  and  for  the  condition  of  time  dependence  e^      ,   the 
force  equation  takes  the   form: 

p(r,e,z)  »  ju>  fo0(r,e,z) 
The   average  pressure   amplitude   on  the   active   area  of  the  h^n  cylinder   is 
defined  as  p^. 

v/oc  (2h-l)bfB 

j"ifi_       \  \  0(a,e,z)   de  dz 


4<XJ9 


-ex       (2h-l)b-yS 


(2.15) 


The  averaging  prooess  is  carried  out  and  the  expression  for  the 
average  pressure  is,  upon  substitution  for  An: 


00 


n-0 


sin  ^n**.) 
;n      2 
(n*} 


H  (2)>a) 
sin  (2  Vfe  )  n   7 

*Z  -(2) 


(Zirff*    )       /*.  JH^/a)  -  H^Oa)] 


-oo 


oo 
x^>  V.  e^ie-J4^^-  ^'h)b  dj° 
i-1 

(2,16) 

The  average  radiation  impedanoe  per  unit  area  of  the  active  area 
of  the  h™  cylinder  will  he  defined  as  <z>  ,  •        v 


<Z>h  = 


Sl 


\   j^h 

h  e 


Further  define 


_1_ 
k 


<r  ■    —  f- 


<r  =  ^1  -   (A?  )2'        (xj>  )2  <  x 

r:    j-   l|(  *Jf  )2  -  i  (Vf  )2>i 


(2.17) 


(2.18) 


Then  the  radiation  impedance  per  unit   area  of  the   active  area  of  the 
h       oylinder  may  be  written: 


oo 


<z>h^4ebc<^ 


n=0 


sin   (na) 


■n 


?  (2) 

sin   (2Tr$fi  )     H?     (fka) 


(n*)2  \  &*W  fr^Vka)  ^((ncaj 


m 


V, 


V, 


JGPi-rh)    -j4>jb  (i-h)b 


r 


i=0 


(2.19) 
The  radiation  impedance   of  an  array  of  finite  cylinders  which  pulsate 
uniformly  with  arbitrary  amplitude   and  phase   over   a  sector  of  each  cylin- 
der has  been  determined.      Evaluation  of  equation   (2.19)   involves  numerical 
integration  of  an   infinite   integral.      The   integrand  of   (2.19)  must  be 
carefully  examined  in  the  region  of  jf-*-  — .     A  further  handicap  to 
evaluation  of   (2.19)   is  the  fact  that   second  order  Hankel  functions   of 
complex  arguments  must  be  evaluated.      These   qualities   of   (2.19)  are 
especially  undesirable   if  tho  radiation  impedance   is  to  be  evaluated  by 
digital  computer  techniques.     Equation   (2.19)  will  now  be  reduced  to  a 
form  suitable  for  numerical   integration. 

The   following  relation  is    stated  and  used  in  deriving  a  new  form 
of  equation   (2.19).4 

Kv(jx)  «  -  Jfre-tiV  (v-1)  Hv(2)(x) 

1 


oo 


2    h=o8fbG  ^jr 


n-0 


•n 


sin    (n<t) 
(n«)2 


(2) 


sin2    (2  7r,f  ye  )       ?W        ((Mca)  x 

\  o      r     (2)  (2) 

{WJP  )V  [Vl  (^ka)-IIn.1((Mca)J 


m 


i=l 


vi  cj(^i-yh) 


cos  [4fry  b(i-h)"|        djf 
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%I,W.  McLachlan,    Bessel  Functions   for  Engineers,   Oxford  at  the  Clarendon 
Press,    1955,    p.    204. 


oo 


■/jqe°  %£ 


oo 


n=0 


sin 


InJ 


(n    )' 


sin2(27)-j'/3  )(-j)  Kn(.i<Tka)  x 

(2^/3  )2  <r  [Kn^1  (jO-ka^Kn.iC jjT-ka)] 


m 


xyj^eJ'(^^h)0osr4nb(i.h)l     df 


i=l 


The   integrals  will  now  be  considered  separately.      Let 


(2o20) 


1_ 


II  = 


sin2(2  Tr^^g)       42)((Tka) 


(2)  (2) 

(2^J^)2  <T  [l^O-ka)   -  Hn-1(<Tka)] 


Let  XJ°  s   sin  fc 

then  when     j»  =  0   ;       i  =  0 


-  cos  [4^  (i-K)bJ  d£ 


(2.21) 


^  *      ^1  -    (  *.£  )2'    =  cos  J 
d£  a       oos  £ 


A. 


I 


Therefore: 


** 


sin   (k/3  sin  ?  ) 
II  =      V    C *- 


IIn(2)(ka  cos  j  ) 


(k/3  sin/  )  ^^(ka  cos  /  )  -  11^  (ka  cos/  ) 


cos  [  2kb(i-h)sinf]df 


A 


(2.22) 
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Consider : 


in'(2frf/3)    (-.1)  Kn(,i<rka) 


;os    [4TTj(i-h)bl    d^ 


Let  A.f  a     csc  fc 

Then  when  J°  =       i  ;       f a  "^ 

T  =     co  ;     /  =  0 


*~  =     -j 


Jujr  )2  -T1  =  -j  cot | 


d^f  =     -  csc£     cot£     — — 


Therefore : 


i„  = 


Kn(ka  cot/  ) 


(2.23) 


sin    (k  jB  cso  (  ) 

(k/>  )     cso/  Kny!  (ka  cot  /  )  -^    Kn-1(ka  cot  j-'  ) 


cos  r2kb(i-h)csc^]d| 


7L 


(2,24) 

The  expressions  (2.22)  and  (2.24)  are  substituted  in  (2.20)  and  the 
resulting  equation  for  radiation  impedance  per  unit  area  z  ^ takes  the 
more  tractable  form: 


<z>  h  a  j8(O0c   j§ 


/- 


ocB 


oo 


n=0 


(2), 


•    2/        \     \       •    2/i  „    •     k\     H   ^'(ka  cos  £  ) 

sm6(n«:  )     1    sinc(tysin  f)       n        v  i  ' 

(2)  (^ 

Hj^V^ka  cosp-Hjj^Cka  cos£) 


(k/Jsin  1)  ^(2) 
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i=l 
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■/■  J8foc 


oo 


n=0 


L.ej(^i-fh)  cos    j^kb(i-h)sin[]dj 
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>     -i- 


,j(Vi  -^h)  cos    [2kb(i-h)csc^J      d^ 


i=l 


(2.25) 
The  evaluation  of  <z>  h  nas  been  reducod  to  that   of  evaluating  finite 
integrals   and  functions   of  positive  real  variables   for  which  mathematical 
tables   are  readily  available.     A  further  manipulation  of    (2.25)  will  be 
made   in  order  to   separate  the  real  and  imaginary  components   of  < z>  ^» 
Let 

Hn(2)(ka  cos  |  ) 


Bn(ka  cos/  )  e"J^n(ka  cos£  )  - 


II. 


(2) 


,2)  (2)       , 

^^(ka  cos  f  )  -  Hn_]_(ka  cos|  ) 

(2.26) 
Then  the  expression  for  < z>  ^  transforms  into  the  following  expression: 
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The  radiation  impedance  per  unit  area  of  an  array  of  m  finite 
cylinders  -which  pulsate  uniformly  over  a  given  sector  of  each  cylinder 
with  phase  and  amplitude  arbitrary  but  constant  for  a  given  cylinder  is 
expressed  by  equation  (2.27).   There  are  no  restrictions  stated  on  the 
wavelength  of  the  radiated  sound,  on  the  cylinder  length  or  radius,  or 
on  the  cylinder  spacing.  Although  equation  (2.27)  is  for  a  system  with 
uniform  velocity  and  phase  across  each  radiating  element,  the  method 
of  solution  is  not  limited  to  such  a  velocity  distribution.  Any 
velocity  distribution  in  azimuth  which  can  be  expressed  in  a  Fourier 
Series  is  possible  of  solution  using  such  technique.   Further  any 
velocity  distribution  in  the  z  coordinate  for  which  a  Fourier  Transform 
exists  is  also  capable  of  yielding  a  solution.   The  distribution  in  the 
z  coordinate  need  not  be  the  same  on  each  cylinder  for  a  solution  could 
be  constructed  piece  by  piece  for  such  a  system. 
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3»  Special  Gases 

The  validity*  of  equation  (2.27)  will  be  considered  by  investigating 
two  special  cases  of  the  finite  array  for  which  solutions  already  exist* 
The  first  array  to  be  considered  is  the  single  uniformly  pulsating 
cylinder  in  an  infinitely  long,  rigid  baffle 9     Let  the  array  be  con- 
structed as  follows: 

1,  There  is  one  active  cylinder  of  length  -  2  f>    and  radius  =  a, 

2,  The  cylinder  is  mounted  coaxially  in  an  infinitoly  long,  rigid, 
cylindrical  baffle  of  radius  =  a. 

3,  The  cylinder  pulsates  uniformly  with  radial  velocity  - 


V0  e- 


,JVo 

The  Fourier  Series  equivalent  of  the  third  condition  as  stated  above 
is  expressable  in  one  term.  That  is,  on  the  active  cylinder  the  radial 
velocity  may  be  expressed  as  one  term  of  the  Fourier  Series  (2,01), 

vr  3  VQ  eJ^o 

With  reference  to  equation   (2.27),   the  above   statement   or  condition 
reduces  the   summation  with  respect  to  n  to  one  term,   namely  n  =  0«     The 
condition  of  one  cylinder  reduces  the   summation  with  respect  to  i  to  one 
term,   namely  i  =  0o     Let  the   radiation  impedance  per  unit   area  of  the 
uniformly  pulsating  cylinder  be    <  z>  .      Then  from  equation   (2.27), 

z       =  8foC    A    \       g£Lik/S  sin  -,  )   Bo(ka  cos  /  )   sin  ^0    (ka  cos      )  d  W 
4    J  (kasin?  )2  ^         •  A 


' 


2 


/    ^Co°  £  \       Sin2(VSlnP   B0(ka  cos  £  )  cos   TQ    (ka  cos  t  )  d  £     .         + 

J         (ysinp2 
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,  £    \        sin2(kBcsci:     )     X0(ka  cot  p) 


0 

(3.01) 
Equation  (3.01)  is  written  with  substitutions  according  to   (2.26): 


<z>  =  j4/onc  A  \    sin2(¥sin(    }     Sa(2^a    003i)  d^     ^ 


(k/3  sin  |  )  l^^^ka  cos*  ) 


A   I       S        (Mo so/)  K0(ka  cot  i  ) 

+  J'4<°oC  %  J  (k/3)2csc|  K1(ka  cot  t)  *  ' 


0 


Equation   (3.02)  compares  with  Robey's  equation   (18)     for  the 


(3.02) 


radiation  impedance  of  a  single  uniformly  pulsating  cylinder  in  an 
infinitely  long,  rigid  cylindrical  baffle  if  the  variables  of  integration 
are  changed  and  it  is  allowed  that  Robey's  equation  contains  an  error  in 
sign.   (This  is  demonstrated  in  Appendix  II). 

It  will  now  be  shown  by  uso  of  the  Dirac  Delta  Function  that 
equation  (3.02)  will  transform  into  the  well  known  expression  for  the 
radiation  impedance  of  an  infinitely  long,  uniformly  pulsating  cylinder. 
Equation  (3.02)  is  written  in  the  form  of  equation  (2.19)  by  trans- 
formation of  variables; 


°D.  H.  Robey,  On  the  Radiation  Impedance  of  Finite  Cylinders,  J.  Acoust. 


Soc.  Am.,  27,  pp.  706*711,  1955 
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sin2(27r/3f )     H0(2Vka)         .  ^ 
<  z>     -   ipc         \         /    *         df 

°     J    ^H,        (CTka) 


Consider  the   expression 

lla  1        sln2(2^Jr> 

2»73*oo  2  7/0  (f  )2 


(3.03) 


(3.04) 


Expression   (3.04)  has  the   following   qualities: 

1.  It   is    zero  for  all      j°    except     y   -   0. 

2.  At     Y  =   0,    (3o04)  tends  to   infinity. 

3.  The   integral   of   (3.04)  from     J°    equals  minus   infinity  to 
V  equals  plus   infinity  is   equal  to  unity. 

The   above   is  the  definition  of  a  Dirao    Itelta  Function.     Therefore: 

2, 


L.3"12??':  b(|p 


lim 
2lrp*oo  ^         2V/3(f)2 


(3.05) 


and 


,00 


H   (2)    (fl-ka) 


lim  <z>     =  j*  o      \  <b  (  f      2_ d  f 

2m3-*oo  \  (2)    ,  N  J 

(TK^        (<Tka) 


-00 


HQ(2)(ka) 

lim      <z>    =  o^oc  —~ 

2"p-*oo  Hi  ^(ka) 


(3.07) 
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The  equation  (3.07)  which  is  the  radiation  impedance  per  unit  area 
of  an  infinitely  Ions;,  uniformly  pulsating  cylinder  compares  with  the 
well  known  expression  for  such  radiation  impedance.6  This  concludes  the 
discussion  of  the  validity  of  equation  (2.27). 


L.  Beranek,  Acoustic  Measurements,  John  Wiley  and  Sons,  Inc.,  1949, 


p0  60 
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4,   Beam  Steering 

The  application  of  equation  (2.27)  tojspecific  array  will  now  be 
considered  in  order  to  study  the  effect  of  electrical  bean  steering  on  the 
radiation  impedance  of  the  array.  The  array  is  to  be  constructed  as  follows: 

1.  There  aro  nine  active  cylinders  with  kp  =  jg-Tf,  ka  =^  and 
kb  s  0.61T' 

2.  The  cylinders  are  mounted  coaxially  in  an  infinitely  long, 
rigid  cylindrical  baffle  with  ka  =  If 

3.  The  cylinders  pulsate  uniformly  over  their  surface  (symmetry 
in  azimuth). 

It  is  proposed  to  investigate  the  effect  of  electrical  beam  steering  on 
the  radiation  impedance  of  the  individual  cylinders.   Electrical  beam 
steering  will  be  accomplished  by  phase  leading  each  cylinder  a  uniform 
number  of  degrees  (  iii>  )  with  respect  to  its  position  in  the  array. 

Equation  (2,27)  is  written  to  conform  to  the  above  boundary 
conditionso 


P2 


<z> 


h  =  j 


<°o° 


-L  (2) 

sin2(i-irsin()^  e^i-h^coSri.27y(i-h)sin>lHo   »QQBf) 


9 


-/  J  \     sin2(^cscO  y  ej(L-*)V      r  VOrcot|)        c 

\  /ir^2_    .    ^— -  cos     1,2^  (i-h)cscf  hT7~ — ^  J  >  i 


(4.01) 
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The  result  of  the  evaluation  of  (4,01 )  by  numerical  integration 

using  the  trapezoidal  rule  and  ^  £   =  — —  is  plotted  in  figures  2  and  3 

■■>     40 

for  cylinders  number  one,  three,  and  five* 

For  such  an  array,  the  electrical  beam  steering  has  but  a  moderate 
effect  on  the  radiation  resistance  of  any  element.  Electrical  beam 
steering  affects  the  radiation  resistance  of  the  end  elements  of  the  array 
to  a  greater  degree  than  it  affects  that  of  the  center  elements. 

Electrical  beam  steering  has  a  more  drastic  effect  on  the  radiation 
reactance  of  such  an  array.  Radiation  reactance  in  this  array  varies  as 
much  as  ten  to  one  for  the  center  element  of  the  array. 
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APPENDIX  1 


The  Helmholtz  Wave  Equation  is   restated  in  cylindrical  coordinates: 

-^f   -f±^-fh^+    ¥*S   -/^=0  (A1.01) 

dr^  r      or         r6     <*  e6  a  z6 

The  solution  to  the  above  equation  is  assumed  to  exist  in  a  variables- 
separable  form: 

jrf(r,e,z)  =  R(r)9(e)Z(z)  (A1.02) 

(Al,02)  is  substituted  into  (Al.Ol)  and  the  resulting  equation  is: 


d2!*   ,  1   dR   ,  1    d2£     .    1   d2Z 


2 


^-  —  ^b  — f^—  ^Hr  ^k  "  °  (A1-03) 

R       dr2         Rr     dr         fir2     d©2         Z       dz2 

1       d2Z 
It  can  now  be   seen  that  —      "  -   is   equal  to  a  quantity  which  is   independent 

Z   dz2 
of  z;   therefore  the  following  is  true: 

-L  -Sff  =  -  (2*JT  )2  (A1.04) 

Z   dz2 

Where   f  is  an  eigenvalue  to  be  determined  by  the  boundary  conditions 

assigned  to  the  array.  The  solution  for  Z(z)  may  now  be  written. 

Z(z)  =  a.    eJ2lr^z  /bY    e   -J2'T^   z  (A1.05) 

1  p 

1        d2Q 
Continuing  with  the   solution,    it  can  be   seen  that  — —  m     is  equal 

©       de^ 
to  a  quantity  that   is   independent   of  e;   therefore: 
I2©  - 


-1-  -dfS  =  _  n2 

e  d© 


(A1.06) 

Where  n  is  an  eigenvalue  to  be  determined  by  the  boundary  conditions 
assigned  to  the  array.   It  will  be  found  more  convenient  to  write  the 
solution  for  ©(e)  as  follows: 

©(e)  ~   cncos(ne)  ■/  d^   sin(ne)  (A1.07) 


23 


The  assignment  of  arbitrary  eigenvalues  to  the  solutions  for  Z(z)  and  ©(©) 
determine  the  form  of  solution  -which  Rn,f  (r)  must  assume, 

d\y  jl  ***y   f[*2  _  (2**  )2  ]  -  sL  }  Vy  =  °      (Alo08) 

dr     r   dr  r 

In  the  following  discussion,   the   quantity  k2  -    (2  7>^   )2  will  be  defined  as 
equal  to./''-2.     The   quantity/^  must  reduce  to  k  as    (2  ^^f  )     approaches   zero; 
therefore : 


y^  =  ^k2  -    (2^^  )2  (27rr)2   <    k2 

^=  -j{(2^)2-k2'      (277J)2  >  k2  (A1.09) 

To  satisfy  the  Sommerfeld  Condition  of  Radiation,  the  solution  for  Rn,f  (r) 
must  be : 

^n,y  (r)  =  Hn(2)  j^r)  (A1.10) 

The  expression  for  Rn  «  (r)  as  r  approaches  infinity  is: 


Rn,y  (r)  =  y 


2_  Q-jy  r  -1h2^i)         (Al.ll) 


y/^  r 
This  expression  coupled  with  the  time  dependence  eJ*4^  represents  an 

outgoing  wave. 

The  wave  equation  has  been  solved  consistent  with  the  Sommerfeld 

Condition  of  Radiation.   It  is  now  possible  to  state  a  eigenfunction  which 

is  a  solution  to  the  wave  equation;  and  if  the  boundary  conditions  imposed 

on  the  array  are  consistent,  a  sum  of  the  eigenfunctions  should  satisfy 

the  additional  boundary  conditions. 

The  eigenfunction  is  stated: 

^  (,...,)  .  AM  ££j.*W  *  H^ r)  (A1.12) 
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APPENDIX  II 
To  substantiate  the   statement  that  equation   (18)  of  Robey's  paper   is   in- 
correct,  the  development   of  equation   (18)  from  equation   (16)  -will  be   inves- 
tigated ap  well   as  the  validity  of  equation   (16).      Further   it  will  be   shown 
that  equations    (19)   and    (20)   of  Robey's   paper  are  correct. 

.-,00 

■qq  "   J*W  Ki^y(JL) -    dT  (16) 

0 
Where  K  =    (k2-J  2)2 
Rewrite  equation   (16)  as   follows: 

-/oo 

<2qq>=.jii^L       /         O&teL     _ik|!iilM     d/(A2„01) 

^  4'raL  J  (ff  KHl(2)(Ka)         J 

-OO 

Where    <zqq>    -"the  average  radiation  impedance  per  unit   area 
Consider  the  following  expression: 

lim        ^n2(^L)  (A2.02) 

L-»oo  ^L(J°)2 

The  following  facts  are  true: 

1.  (A2.02)  approaches   infinity  as  j°     approaches   0, 

2.  (A2.02)   is   zero  for  all  other  f  . 

3.  The   integral  of    (A2.02)  from  f  =  -00  to   J3   =  -/co  is  equal  to  unity. 
Therefore    (A2.02)   is   a  Dirac   Delta  Function,    and   (A2.01)  may  be  written 

as  follows : 


H0(2)(Ka) 
L->oo  X    qq'  *"*"      /  VJ    '  K  II1(2)(Ka) 


lim      <z     >a-jjf0W  b(f)   °  *T  (A2.03) 


-00 


26 


Thus  the  radiation  impedance  per  unit  area  of  the  uniformly 
pulsating  cylinder  as  derived  by  Robey  reduces  to  the  following  as  the 
length  of  the  cylinder  approaches  infinity: 

H0(2)  (ka) 
1dm  <z   q>  a  -jfto  -2— (A2.04) 

L"°°  KX{Z)    (ka) 

This  does  not  agree  with  the  well  known  result  for  an  infinite 

H  (2)(ka) 
cylinder.  Further  the  angle  associated  with  mQ/r,\   -"   ^s  always  in  the 

V   (ka) 
fourth  quadrant;  thus  the  radiation  resistance  given  by  (A2»04)  would 

always  be  negative. 

It  can  be  assumed  that  Robey1  s  equation  (16)  is  correct  except  for 

sign,  and  it  will  be  shown  that  Robey* s  equation  (18)  has  an  error  in 

sign  but  that  equations  (19)  and  (20)  are  oorrect.  Equation  (16)  is 

restated  with  the  correct  sign: 

■„  -  w  /   a4S*  ^r^  *r         (--a) 


0 


W^        ml    W   (Ka) 


k 


sm 


2   (j>L)  H0(2)(    fk2-?2     a) 

(JT)   2  fk2-?2'   II1(2)(|k2-^     a) 


zqq  =  38a^  I    — -rrr-   r^-r^  ,.  ,^t  ■  ,»    ^    ,      df 


sm 


2(J°L)  (-j)£0(   ^J2-k2'     a) 


dr 


k 


{ff  (-j)ff2^    K^jT^F^ 

(A2.07) 
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To  reduce  (A2.07)  to  the  form  of  Robey's  equation  (18)  consider  the 
following  integrals: 

.k 

II  = 

0 
Let      "fa  k  sin  -&■ 
then  when    1  =  0,    ~&  =  0 


si*2(y  L)  "o(2)(     ^2-y2       a)  ^ 

<n2     iXFV5VT7~rra) J 


7  n  k,    ty-a  ^ 

9 

and 

]|k2  -f2'  -  k  cos^ 

d]f  a  k  cos-z^  &zy 

Therefore 

ft 
C2 

h* 

1           \           sin2(kL  sin-^)       H0(2)(ka  cos  ■&) 

k2        J                 (sin^)2             H/^dcaooatf) 
0 

Consider 

oo 

X2= 

k 

sin2CfL)                K0(i 

J"2  -k2     '       a) 

(J)2                  2"k2 

K1(^2-k2'     a) 

d^ 


Let  $  =  k  cosh/' 

then  when     ^  =  k,   f  »  0 
*  =  go,  fa  oo 
and    jj>2-k2     =  k  sinh<^ 

&f  =  k  sinh^/  d<f 


(A2.07) 


(A2.08) 


dj  (A2.09) 
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Therefore 


I     -  JL        \         sin2(kL  cosh^  )     KQ    (ka   sinh  f>   )  ^  (A2.10) 

Kx    (ka  sinh^   ) 


(cosh^)2 


Equation   (A2.05)  may  now  be  written  as 

IT 


2  (2) 

ppCQ  \        sin    (kL  sin^)  IIQ        (ka  cos^)  , 

k? 


99  -  J         p  \  2         (jo 

J  (sinx*)         H^    ;(ka  cos^) 


,      8a.^0u>         \         sin2(kL  cosh^)  Kn(ka  sinhf  ) 
-/j_£o_ ^       jv   _2 dV,  (A2.ll) 

k  o  (cosh/'  )  K1(ka  sinh/-'') 

Thus  from   (A2.ll),    it  can  be   seen  that  the  terms   in  Robey's  equation 
(18)   should  have   identical   signs;    and  from  the   initial   argument   about 
equation   (16), •  the  terms   should  be   positive. 

Robey's  expressions  for^and^  as  they  appear  in  equations    (19)   and 

(20)  will  now  be  derived  from   (A2.ll).     The  &  torm  is   self  explanatory  and 

correct.     Examine  the   expression 

Hn(2)(kacos^)        m     Koiilix)  where  x  =  ka  cos  ^ 

Il/^ka  costT  )  r^Z)    (x) 

J0(x)  -  jYo(x) 
Jx(x)  -  jY-^x) 

s     Jnto  -  .1Yn (*) £li*L£j!li*L 

Jx(x)  /Y1(x) 
-     Jq(x)Ji(x)   •/Y0(x)Y1(x)  J0 (x)Yj(x)  -  J]»Y0(x) 

J2(x)   /y*(x)  J1(x)-/Y1(x) 
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J0(ka  cos  ^)Ji(ka  cos  z.*)  /-Y0(ka  oost^)  Yi    (ka  oos^) 

2  2 

J    (ka  oosv)   -/Y  (ka  oos  ty) 


-  J 


ka  cos  i*"  .2,,  ,      /    2,  v 

J    (ka  cos^j  ■/•Y .(ka  cos^) 

Thus   it  can  be   seen  that  Robey's  coefficients  /c' ,  $  ,   and  <$  in 
equations    (19)   and   (20)   are  correct   as   stated  but  are  not  consistent  with 
equation   (18 )o 

The   above   arguments   are   believed  to  be    sufficient  to   show  ilobey's 
equation   (18)   is  not  correct   as   stated,   but  would  be  correct   in  the  form 
of  equation   (A2.1l)»     Equation   (A2.ll),    if  stated  in  terms   of  average 
radiation  impedance   and  with  a  change   of  variable,    is  equivalent  to 
equation   (3.02). 


7 

Go   N.  Watson,   A  Treatise   on  the  Theory  of  Dessel   Functions,   The  Macl.lillan 

Company,    1945,    p   77 
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(Figure  ±  ) 
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